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This is a short review of recent developments in calculation of higher order corrections to various two-point core¬ 
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1. Introduction 

Precise determinations of parameters of the Standard 
Model (SM) and precise predictions for observables 
measured at present and future experiments are critical 
in testing the SM and may hint towards physics beyond 
the SM. Increasingly precise results from high energy 
experiments at LEP, LHC or a future electron-positron 
collider have been obtained during the past years or are 
expected for the coming decade, with production and 
decay rates or masses of gauge or Higgs bosons or of the 
top quark as characteristic examples. These are comple¬ 
mented by measurements at low energies, which lead to 
precise values of the strong coupling from r lepton de¬ 
cay or the masses of strange, charm and bottom quarks. 

To extract the fundamental parameters of the theory 
and relate the large number of experimental results, the 
knowledge of higher order perturbative corrections is 
crucial. Significant advances in this direction have been 
made during the past years, in particular in the frame¬ 
work of the collaborative research center “Computa¬ 
tional Particle Physics” (SFB/TR-9). 

Within perturbation theory quantum-theoretical am¬ 
plitudes are described by Feynman Integrals (FI’s). Im¬ 
proved precision, which is required both for strong and 
electroweak interactions, necessarily leads to a signifi¬ 
cant increase of the complexity of the calculations. This 
applies to the number of FI’s, their increasing complex¬ 
ity and, consequently, to the effort required for their 


evaluation. 

The complexity of a FI can be roughly measured by 
the sum of two quantities; (i) the number of (indepen¬ 
dent) external momenta and (ii) the number of so-called 
“loops”, that is the number of integrations with respect 
to internal momenta which should be performed. In ad¬ 
dition, the pattern of masses of (virtual) particles ap¬ 
pearing inside of the FI also presents an important fea¬ 
ture characterizing the complexity of the integral. 

For instance, one-loop diagrams can be calculated 
in analytic form for arbitrary masses and external mo¬ 
menta, diagrams with three or more loops, however, 
only for configurations involving just one mass or en¬ 
ergy scale (one-scale diagrams). Problems involving 
several scales, in particular those with pronounced scale 
hierarchy, can be treated approximately using asymp¬ 
totic expansions (Hard Mass Expansion and/or Large 
Momentum Expansion, for example), where each of the 
diagrams is expressed through a nested sum of one-scale 
diagrams HI. 

The present review will deal with a special class of 
one-scale FI’s, namely massless propagators, that is in¬ 
tegrals depending on only one external momentum, q, 
and with vanishing internal masses. In what follows we 
will customarily refer to massless propagator-type FI’s 
as p-integrals. 
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2. Massless Propagators and Physics 


2.1. RG functions and IR reduction 


The method of the renormalization group (RG) 
is of vital importance in modern quantum field theory. 
It is enough to recall that the famous idea of asymptotic 
freedom is based on the RG concept of the running cou¬ 
pling constant. 

The RG functions — jS-functions and various anoma¬ 
lous dimensions — serve as coefficients in the RG equa¬ 
tions. They can be conveniently expressed in terms of 
p-integrals (see below) within the framework of Dimen¬ 
sional Regularization SO and Minimal Subtraction 
(MS) schemes IH. The naturalness and convenience of 
the MS-scheme for RG calculations comes from the fol¬ 
lowing statement 0: 

Theorem 1. Any UV counterterm for any FI integral 
and, consequently, any RG function in an arbitrary min¬ 
imally renormalized model is a polynomial in momenta 
and masses. 

This observation was elaborated by A. Vladimirov QO) 
to simplify considerably the calculation of the RG func¬ 
tions. The method was further developed and named 
Infrared Rearrangement (IRR) in m. It essentially 
amounts to an appropriate transformation of the IR 
structure of FI’s by setting zero some external momenta 
and masses (in some cases after the differentiation is 
performed with respect to the latter). As a result the 
calculation of UV counterterms is reduced to that of p- 
integrals. The method of IRR was ultimately refined 
and freed from unessential complications by inventing 
the so-called R*-operation IfT^ . The main use of the 
R* -operation is in the proof of the following statement 

Cl: 

Theorem 2. Any (L+l)-loop UV counterterm for any 
Feynman integral may be expressed in terms of pole and 
finite parts of some appropriately constructed L-loop p- 
integrals. 

Theorem 2 is a key tool for multiloop RG calculations 
as it reduces the general task of evaluation of (L-i-1)- 
loop UV counterterms to a well-defined and clearly 
posed purely mathematical problem; the calculation of 
L-loop p-integrals. In the following we shall refer to 
the latter as the L-loop Problem. A short account of the 
current status of the Problem can be found in Section[3 


2.2. Two-point correlators 

Within perturbation theory, every two-point correla¬ 
tor 

W'j- = Jdve''?-' <0|r [j2(x) ;i(0)] |0), (1) 

with ji and jz being in general elementary fields or 
(local) composite operators, is expressed within PT in 
terms of p-integrals provided the momentum transfer q 
is considered as large with respect to all relevant masses 
and, thus, the elementary field propagators contributing 
to YV'ti be effectively considered as massless. 

An important class of two-point correlators is repre¬ 
sented by the case of ji and 72 being quark currents of 
the form: 


ii = ipTfi, j2 = j\. 


In particular, the total cross-section of e^e^ annihila¬ 
tion into hadrons, the (inclusive) Higgs decay rate into 
hadrons, the semihadronic decay rate of the t lepton 
coupling are all expressible in terms of absorptive parts 
of the quark current correlator Q with r chosen as yy,, 
1 and (1 - 75 ) 7 ^ respectively. 

Clearly, one could compute a (LH-l)-loop two-point 
correlator by computing the corresponding set of (L-i-1)- 
loop p-integrals. But one can do much better if the final 
aim is the absorptive part of the correlator. 

Indeed, let F be a particular (L-i-l)-loop Feynman 
diagram contributing to the perturbative expansion of 
a massless correlator. The renormalized version of 
the corresponding Feynman integral can be generically 
written a^ 


R<r>((2') = <r)(e2)- 


ZyZy(r/7}(Q^) + ... ( 2 ) 


Here Zy is the UV Z-factor corresponding to a IPI sub¬ 
graph 7 of F and dots stand for contributions with two 
and more UV subtractions. The finiteness of the left 
part of eq. Q means that the pole part in e = (4 - D)/2 
of (r)( 2 ^) is completely fixed by poles in e which ap¬ 
pear in UV subtractions (the boxed terms in Q). On the 
other hand, the UV subtractions could, obviously, con¬ 
tain L'-loop Z-factors with L' < L -H 1 and the reduced 
p-integrals like (F ly){Q^) with the loop number not ex¬ 
ceeding L! Applying Theorem 2 we arrive at the con¬ 
clusion that the pole part of (F)( 2 ^) (and, consequently. 


'without essential loss of generality we assume that {T)(as, Q^) 
is a scalar integral depending on the external momentum Q via its 
square, Q- = QvQ'■ In addition, we set the renormalization scale 
parameter^ = 1. 
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its absorptive part) is completely expressed via L-loop 
p-integrals only. 

We want also to stress that by high-energy limit we 
understand not only the case when all masses can be 
neglected but also the possibility to take into account 
mass effects by exploiting a small mass expansion. As 
a suitable example one could mention the calculation 
of the power suppressed (of order ni^ls^ and so 

on) corrections for the correlators of (axial)vector quark 
currents in higher orders of pQCD lUMIl. 

2.3. OPE and DIS 

The theoretically cleanest description of (inclusive) 
Deep Inelastic Scattering (DIS) can be achieved within 
Operator Product Expansion (OPE) of two composite 
operators (for a recent review see, e.g IfTSl ). Here the 
main objects to compute are the so-called Coefficient 
Eunctions (CE) which can be always computed via p- 
integrals with the help of the so-called method of pro¬ 
jectors imiioi. It is important to stress that within the 
method of projectors one needs no IRR: L-Ioop correc¬ 
tions to a CF can be expressed directly in terms of L- 
loop p-integrals. 

A good example of an early multiloop OPE calcula¬ 
tion is the one of the corrections to the Bjorken sum 
rule for polarized electroproduction and to the Gross- 
Llewellyn Smith sum rule 11211 . We will discuss later 
our calculations of the next, order contributions to 
the Bjorken sum rule. 

3. Calculational Methods 

A signihcant number of higher order calculations are 
usually performed according to the following “stan¬ 
dard” scenario. First, the Feynman amplitudes are re¬ 
duced to a limited set of so-called master integrals 
(Mi’s). This step is based on recursion algorithms 
obtained by using Integration-by-Parts (IBP) identities 
(see, e.g. recent books and reviews lEHlSl and refer¬ 
ences therein). 

An important feature of the standard scenario is that 
the resulting set of master integrals should be computed 
only once and forever due to the well-established prop¬ 
erty of universality: for every given class of Feynman 
amplitudes characterized by the number of loops and 
the pattern of external momenta and masses the corre¬ 
sponding set of master integrals is universal in the fol¬ 
lowing sense: every (even extremely complicated) am¬ 
plitude from the class can be expressed in terms of one 
and the same (finite! Il26ll27l l set of master integrals . 

Thus, the task of evaluation of p-integrals at L-loops 
(L-loop Problem) is naturally decomposed in two: (A) 


reduction of a generic L-loop p-integral to masters and 
(B) evaluation of the latter. 

Both A and B Problems were solved at two- and 
three-loop level long ago MM- The four-loop prob¬ 
lem has been under active investigation in our group 
since the beginning of the current century. We will de¬ 
scribe the current status in the next two Sections. 

3.1. A: Reduction via 1/D expansion 

The standard (Laporta) Il29l |30| approach to solve 
IBP relations implies a step-by-step linear reduction of 
more complicated integrals to less complicated and fi¬ 
nally to irreducible (master) integrals. Unfortunately 
this conceptually simple method could not be used in 
our case because of the extremely large amount (10^ - 
10*) of 4-loop integrals appearing after IRR for a typical 
5-loop problem. As the result we used a more sophis¬ 
ticated and laborious, but less demanding for computer 
resources method based on large D expansion ll3Tll of 
the formal solutions of the IBP relations ||32l. 

So assume that we need to perform the reduction, that 
is to calculate the coefficients in front of master inte¬ 
grals. The coefficients depend on indices of the original 
integral and fulfill the IBP identifies according to these 
indices. If we construct some “convenient” solutions of 
the identities, then the coefficients we need can be ob¬ 
tained as their linear combinations with proper bound¬ 
ary conditions. 

As shown in 0^ such “convenient” solutions can be 
constructed in the form of the integrals of polynomial 
raised to some degree (linear in dimension D). Unfor¬ 
tunately, although these integrals are simpler than the 
original FI, they are still too complicated for direct eval¬ 
uation. From the other side, we are interested in linear 
combinations of these integrals which are rational over 
D (because in principle they can be calculated by the 
step-by-step reduction). So one can expand these inte¬ 
grals in the 1 /D limit (resulting in integrals of Gaussian 
type in]), calculate sufficiently many terms and finally 
reconstruct the exact D dependence. 

Note, that we need not calculate each of the 10^-10* 
integrals involved in the specific problem individually 
(as it is necessary in the standard reduction). Instead 
we can calculate 1/D expansion coefficients of the to¬ 
tal expression we are interested in, thus saving a lot of 
computer resources. 

The construction of the large-D limit requires in gen¬ 
eral huge storage resources, which naturally constrains 
the structure of the input p-integrals: they should better 
not contain any extra parameters like color coefficients, 
nf, the number of light quark flavours contributing to 
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Figure 1: Non-trivial three-loop master p-integrals. 

internal fermion loops, and so on. As a result we are 
forced to use a “slice” approach: that is to set all color 
coefficients to their numerical values and to fix to 
some integer. Combining together different slices one 
can always reconstruct the full «/-dependent structure 
(and, if necessary, even all colour coefficients) 

3.2. B: Master p-integrals and their evaluation 

At one- at two-loop levels the master p-integrals are 
trivial as they can be easily performed in terms of T- 
functions for a generic value of the running space-time 
dimension D (see, e.g. 1331. At three loops there exist 
only two non-trivial master p-integrals (see Fig. 1); their 
values are known since long im. 

The significantly more complicated problem of iden¬ 
tifying and computing all 28 four-loop master p- 
integrals has been solved only recently. We refer the 
interested reader to the original publications OTl [34]- 

|36l. 

3.3. Computer Algebra & FORM 

Higher order calculations dealing with thousands of 
diagrams already at 3-loop level require heavy use of 
computer algebra tools. We are using QGRAF llJTl for 
automatic diagram generation as well as a collection of 
Mathematica and PERL scripts to automatically assign 
topologies and prepare input files for FORM. 

The workhorse for all the complicated calculations 
discussed in the current paper is the computer al¬ 
gebra program FORM li3^ and its parallel versions 
ParFORM |[39l and TFQRM l40l . The program offers 
excellent possibilities for dealing with gigantic data 
streams generated during the reduction procedure. The 
internal specifications allow FORM to deal with expres¬ 
sion which are much larger than the available memory 
(RAM). The only restriction for the size of an expres¬ 
sion is the disk space which nowadays is rather cheap. 
As a consequence, the complexity of a problem solvable 
by FORM is practically restricted only by time. 

The FORM program, BAICER, intended for the 
reduction of complicated four-loop p-integrals imple¬ 
ments the algorithms described in the previous Section. 


With increasing experience and using heuristic criteria 
about the pole structure of the coefficient functions in 
front of master integrals, BAICER has developed into 
an efficient tool which allows to calculate complicated 
four-loop massless propagator integrals, including their 
finite part. It runs routinely on ParEORM and TEORM 
using 8 to 16 cores with a speed-up between 6 and 12. 

4. Scalar Correlator & -» qq) 

The decay width of the Higgs boson into a pair of 
quarks can be written in the form 

m ^ //) - m](p) {s^ Ml, p) (3) 

4 V27r •' 

where p is the normalization scale and 

R^{s) = Imn‘^^(-i - ie)l{2ns) (4) 

is the spectral density of the scalar correlator 

U^^(Q^)^(4nfi J dxe'‘>\OI T[ \0). (5) 

Here - -q^ and - 'Vf'i’f is the scalar current 
for quarks with flavour / and mass mf, coupled to the 
scalar Higgs boson. The 0(a‘l) result for R^ is known 
analytically since long Eniia (early results for orders 
and can be found in li43l and Il44l respectively). 
Eor brevity we put below the final result in numerical 
form: 

R(s,p^ = s) = 1 H-5.6667aj-H [35.94- 1.359 
+a] [164.14 - 25.77 n/ -H 0.259 (6) 

-H a* [39.34 - 220.9 nf + 9.685 nj - 0.0205 nj]. 

We will discuss the application of (j^ for the dominant 
b-quark decay mode of the Higgs boson later in Sec¬ 
tion |7| 

5. Vector Correlator & R(s) 

The ratio 

R{s) = cr{e*e~ —> hadrons)/cr(e^e“ ^ pV~) 

is expressed through the absorptive part of the vector 
correlator 

^pM) = i f dxe'‘i^(0IT[ ]| 0 ) = 

= i-gpvq^ + qpqv)n{-q^), 0) 
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Figure 2: Lowest order non-singlet (a) and singlet (b) diagrams con¬ 
tributing to the polarization operator. 


with the hadronic EM current y™ = Z/ 

Qf being the EM charge of the quark /. The optical 
theorem relates the inclusive cross-section and thus the 
function R{s) to the discontinuity of fl in the complex 
plane 

/?(s) = 127rlmn(-s-i5). (8) 


Eor the vector correlator the terms of order and a\ 
are known since long mu HD. The next, a] order has 
been under investigation for more then 10 years naEi 
HtHSI in Al group. By now it is known in a complete 
form for a generic colour group G i55]|56l. We put 
below only physically relevant result for G = 5 G(3): 


/?(*)= 3^22 ( 9 ) 



- 114 - 3 - 

■ 17"/ + 
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+ a® [n]{ 
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162 108^ 
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+ 36^ ■ 

, 262 / 
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288 48 
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1103 , 275 

— H3 -H — 



+ [«/( 

6131 , 11 2 

5832 432 

4 - 203 
'■ 324 H3 

54 ^'^3 + 

^ft) 

-Hn). 

( 1045381 

1 15552 

593 

432 

71^ 

40655 r 
864 




11 ^2 , 5 ^2 260 ^ 'i 

+ 12 ^ s3 + 6 b3 27 ^5 j 


-2-2 S3 ^ 432 

144939499 _ 49775 _2 _ 5 
20736 384 ^ 


5445 >-2 659- 

■ 8 b3 281 


-I- 



-io) 

745 , 65 , 5 >-2 25 >- 'i 
432 24 b3 + 6 S3 12 SSj 


, /5795 8245 >- 55 /2 , 2825 / 'll! 

+ IW “ 144 « - T fc3 + ^ , 

where = asjn and we have set the normalization 
scale - i; the results for generic values of p can 
be easily recovered with standard RG techniques. Note 
that the two specific quark charge structures in (j^ cor¬ 
respond to the so-called non-singlet (numerically dom¬ 
inant) and the singlet contributions (see Eig. to the 
vector correlator Q. Numerically, 

R{s) = 3^(2j|l -Ha, H-a?(l.986-0.1153n/) 

-H a^ (-6.637 - 1.200 n/-0.00518 n}) (10) 

-Ha^( - 156.608-H 18.7748 -0.797434 

-H 0.0215161 n})| 

- (^ 2/) (1.2395 a^ -H (17.8277 - 0.57489«/)a^) . 

Specifically, for the particular values of «/ = 3,4 and 5 
one obtains (for the terms of order and we have 
explicitly decomposed the coefficient into non-singlet 
and singlet contributions): 


R"2=3(s) = 2 


1-Ha,-H 1.63980^ - 10.2839a2 (11) 


R"/=4(j)= 


-H a 


-106.8798a 

1 -H a, -H 1.5245aj 
-H a^ ( - 11.686 = -11.52 - 0.16527®') 

-H a^ ( - 94.961 = -92.891 - 2.0703®')j , (12) 

■"H=®(i) = y 1 -Ha, -H 1.40902a2 
-H a® ( - 12.80 = -12.767 - 0.037562®') 

( - 80.434 = -79.981 - 0.4531®')]. 


(13) 


Note that for n/ - 3 the singlet contributions vanish in 
every order in a, as the corresponding global coefficient 
(2/ 2/)^ happens to be zero. Implications of this result 
for the determination of a, in electron-positron annihi¬ 
lation and in Z-boson decays are discussed in ||53]|56l. 


As a by-product of the calculation of R(s) the authors 
of ll54l have obtained the five-loop y6-function in pure 
QED, that is a theory with n/ single-charged fermions 
minimally coupled to the photon field. The result reads 
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(the four-loop result is known since long from lISSl l 
= rif 

+ 


[4 4^1 


44 ,1 

3 

-i-4n/A3 -A"^ 

2«/ + y «/ 


1232 


760 2 832 2 

-46«,+ — 


+ aH nj 


+ nf 
+ «/ 


21758 16000 416 1280 

+ - -^^5 


81 

7462 

4157 


27 3 

• 992^3 H- 2720^5 

f856 128 


128^3 


'V 


H-4'3 

243 27 


(14) 


Here the QED coupling constant 

A(ju) = a(jj)li4n) = eijif I{l6n^). 


charge RC, Zg - 


N 

II 

N 

(18) 

Z, = ^(Z3)-', 

(19) 

Z, = Z^^(Z3)-3/2, 

(20) 

Z, = Z“"(Z3)-'/2(Z3Q-', 

(21) 

Z^ = zf^^(Z3)-'/2(Z2)-'. 

(22) 


Within the commonly accepted MS scheme RCs are 
independent of dimensional parameters (masses and 
momenta) and can be represented as follows 

l<t<i ,/ 

Z(h) = 1 + X (23) 

ij 


6. QCD RG-functions 

Our starting point is the QCD Lagrangian with ny 
quark flavours written in terms of renormalized fields, 
coupling constant g and quark masses mp: 

£o = -^Z 3 id^A, - dyA^f 

-^gZl^d.A'i-dyApiA.xAyf 
- (A^ X _ ^(dyA,f (15) 

H- Z^ dyC (dyC) + g Zf ^ &‘c(Axc) 

Uf 

+ ^ f-'(iZ 2 9 + '«/) 

/=! 

with 

{AxB) = ^ /- 

oxfj 

and with bare gluon, quark and ghost fields related to 
the renormalized ones as follows: 

A^ = V^A“^ cl - ^c“. (16) 


where h - g^/(16;7r^) = aj/(47r) and the parameter e 
is related to the continuous space time dimension D via 
D = 4-2e. Given a RC Z{h), the corresponding anoma¬ 
lous dimension is defined as 


yih) = 


A\ogZ{h) 

A/p- 


^ Z„_i n h" 

n=\ 


X(r)« 


n=0 


(24) 


The anomalous dimension of the quark-gluon coupling 
constant h is conventionally referred to as “QCD /?- 
function”; equation pT) leads to the calculationally 
simplest presentation of the function: 

/3(h) = 2r“*-2 75-73. (25) 


The quark mass anomalous dimension, y^, governs 
the evolution of the quark mass, viz. 


2 d v~i 

B = mym{h) = -m y ^ y; h 


i>0 


(26) 


To calculate y^ one needs to hnd the so-called quark 
mass renormalization constant, Z„,, which is defined as 
the ratio of the bare and renormalized quark masses, viz. 


The vertex Renormalization Constants (RCs) 

Zl', T 6 {3g, 4g, ccg, (A<Ag) (17) 


7 1, 

nif Z2 ’ 


(27) 


renormalize 3-gluon, 4-gluon, ghost-ghost-gluon, 
quark-quark-gluon vertex functions respectively. The 
Slavnov-Taylor identities allow one to express all vertex 
RCs in terms of wave function RCs and an independent 


The final formula for y,„ follows from the QCD La¬ 
grangian (15 1 and reads 


Tm = - 72- 


(28) 

























P. A. Baikov et al. /Nuclear Physics B Proceedings Supplement 00 (2015) 1-\16\ 


1 


Thus, to compute the QCD jS-function and the quark 
mass anomalous dimension at five loops one should 
evaluate five separate anomalous dimensions, viz. 

rf®, rl, 73, g/,/,,72- 

By now we have computed all of them except (most dif¬ 
ficult) 73 . The results are presented in the four next sub- 
section^ 

6.1. Five-loop running of the ghost field 
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6.2. Five-loop running of the ghost-ghost-gluon vertex 
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man — gauge fixing condition. The physically relevant ym andfunc¬ 
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Note that the leading renormalon contribution =» a'/' 

vanishes (in any gauge!) due to the Taylor theorem 
which states, in particular, that s 0 in the Landau 
gauge. 


6.3. Five-loop running of the quark field 
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6.4. Five-loop running of the quark mass 
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Note that in four-loop order we exactl}j^ reproduce 
well-known results obtained in lEiiiTa. The boxed 

^This agreement can be also considered as an important check of 
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«/ 

3 

4 

5 

6 

iymfr^ 

198.9 

111.6 

41.8 

-9.8 

Ga 

162.0 

67.1 

-13.7 

-80.0 

^(7«.)4^"^'' Ell 

163.0 

75.2 

12.6 

12.2 

EH 

164.0 

71.6 

-4.8 

-64.6 


7. Phenomenological applications of 


7.1. RG invariant quark mass 
The solution of eq. (j2^ reads: 


mill) 

m(juo) 


cjasiti)) 
c(asOuo)) ’ 





Table 1: The exact results for (yni )4 together with the predictions made 
with the help of the original APAP method and its two somewhat mod¬ 
ified versions. 


terms in ( |6.4[ ) are in full agreement with the results 
derived previously on the basis of the 1/n/ method in 

lEsca. 

In numerical form y,„ reads 

ym - -as - a^s (4.20833 - 0.138889n/) 

-a] (19.5156 - 2.28412n/ - 0.0270062nj) 
(98.9434 - 19.1075n/ 

+ 0.276 163nj + 0.00579322n}) (54) 

-al (559.7069 - 143.6864 + 7.4824 

+ 0.1083n^ - 0.000085359 n)) . 


Let us compare our numerical result for iym )4 

^(ym)4 = 559.71 - 143.6n/ + 7.4824nJ (55) 

+0.1083 0.00008535 n), 

with an old prediction based on the “Asymptotic Fade 
Approximants” (APAP) method m (the boxed term 
below was used as the input) 


45 


(rm)^^^ = 530 - 143 «/ + 6.61 n] 


(56) 


+0.037 


0.00008535 n) 


c{x) = (.r)^‘»|l +dix + id\l2 + di) (58) 

+{d\l6 + dxd2 + d^) + (4/24 + d\d2l2 

+4/2 + did^ + d^ x^ + (9(x^)^, (59) 

di ^-Pxyi) + yx, (60) 

di -4ro/2-y6270/2-^8171/2 + 72/2, ( 61 ) 

d 3 - -4 ro /3 + 2^1 P2 70/3 -h ro /3 ( 62 ) 

+p\y\l^ -^271/3 -^172/3 + 73/3, 
d 4 ^~ 0 \ 70/4 - 34^82 70/4 + 4 70/4 ( 63 ) 

+ y8i >63 70/2 - j 84 70/4 - p\ 71 /4 


+ ^1/3271/2-183 71/4+4 ’ 72/4 
-16272/4-16173/4 + 74/4. 

Here 7 ,- = iym)ilPo, A = /?,7ySo and 

flias) - - X A' = -^ 0 1 Xi A 

i>0 [ )>0 j 

is the QCD /3-function. Unfortunately, the coefficient d 4 
in eq. ( |63| does depend on the yet unknown five-loop 
coefficient fin (up to four loops the /l-function is known 
from iSlHlTl). 

Numerically, the c-function reads: 

c(x) = X^^'^ Cs(x), c(x) = Cc(x), 

nf=3 nf=4 

c(x) = Cbix), c(x) x"^^^ C,(x), 

«/=5 /i/-=6 


However, this good agreement is broken for fixed val¬ 
ues of nj due to severe cancellations between different 
powers of n as one can see from the Table 1. 


all our setup which is completely different from the ones utilized at 
the four-loop calculations. 


with 

c,(x) = 1 +0.8950 x+ 1.3714x2 (64) 

+1.9517 x^ + (15.6982 - 0.11111/34) x\ 
q(x) = 1 + 1.0141 x+ 1.3892x2 (65) 

+ 1.0905 x2 + (9.1104-0.12000yS4)x^ 

Cbix) = 1 + 1.1755 x+ 1.5007x2 


( 66 ) 
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+0.17248 + (2.69277 - 0.13046y64) 

q(x) = 1 + 1.3980 X + 1.7935 x^ (67) 

-0.68343 x^ + (-3.5130 - 0.14286yS4) x'*. 

Eq. @ could be used to define the important concept 
of the RGI mass 

= m(jJi))lc(a,(fiQ)), ( 68 ) 


A remarkable property of the RGI mass is ji and scheme 
independency: in any (mass-independent) scheme 

lim ni{jx) - 

p-^oo 


Due to this property the RGI mass and function c(x) 
are often used in the context of lattice simulations. For 
example, the lattice ALPHAcollaboration uses (68 i to 
find the MS mass of the strange quark at a lower scale, 
say mj( 2 GeV), from the mf^^ mass determined from 
lattice simulations (see, e.g. E]). For example, setting 
= 2GeV) = - 0.1, we arrive at (here the 

formal parameter h = 1 counts loops): 


m,(2GeV) = (fl,(2GeV))5 x 

(l + 0.0895 + 0.0137 + 0.00195 

+ (0.00157 - 0.000011 ^ 64 ) h^). (69) 


get 


^■^(i = = Mh) = 1 +5.667 a, (71) 

+29.147 aj + 41.758 -825.7 

= 1 + 0.2041 + 0.0379 + 0.0020-0.00140, 


where we set as - as/n = 0.0360 (for the Higgs mass 
value Mh = 125 GeV and a,(Mz) = 0.118). The 
decay rate © depends on two phenomenological pa¬ 
rameters asiMfj) and the quark running mass m^. Fet 
us consider, for definiteness, the dominant decay mode 
H —» bb. To avoid the appearance of large logarithms 
of the type the parameter ji should be cho¬ 

sen around Mh- However, the starting value of mj, is 
usually determined at a much smaller scale (typically 
around 5-10 GeV IfSOl ). The evolution of «&(//) from 
a lower scale io ji - Mh is governed by eqs. (57]|63 


which depend on the quark mass anomalous dimension 
■y(as) and the QCD beta function j8(as) (for QCD with 
«/ = 5). In order to match the Ocas'*) accuracy of (7.2 1 
one should know both RG functions j3 and y,„ in the 
five-loop approximation. 


Fet us assume, conservatively, that 0 < < 200. 

The value of = Mh) is obtained with RG running 
from mh(/u = 10 GeV) and, thus, depends on p and 7 „,. 
Using the Mathematica package RunDe(0 IISTlI and eq. 
•iz© we find for the shift from the five-loop term 


In order to have an idea of effects due to the five-loop 
term in (69 1 one should make a guess about ^ 84 . By in¬ 


specting the available four-loop result 


yS(n/ = 3) = 


(70) 


(a, + 1 .777 a] + 4-4711 + 20.990 a* + A 


we conclude that ^84 = 50 - 100 could serve as a nat¬ 
ural estimate of ^ 84 . With this choice we conclude that 
the (apparent) convergence of the above series is quite 
good even at a rather small energy scale of 2 GeV. On 
the other hand, the authors of El cite an estimation 
^84 = -850(!) for the nf - 3 QCD. With such a huge 
value of pA, the five loop term in (69 1 would amount to 


0.01092 and, thus, would significantly exceed the four- 
loop contribution (0.00195). 


dmlCMn) 

mliMn) 


= -1.3 ■ 10-\p4 = 0) 


(72) 


I - 4.3 ■ = 100)1 - 7.3 ■ = 200) 

(73) 


If we set = Mh, then the total effect of OCa*) terms as 
coming from the five-loop running and four-loop con¬ 
tribution to on T{H —> bb) would be around -2%o 
(for pA - 100). This is to be compared to the para¬ 
metric uncertainties coming from the input parameters 
asCMz) = 0.1185(6) ED and «*(;«*) = 4.169(8) GeV 
m which correspond to + l%o and + 4%o respectively. 


8. Deep Inenalsic Scattering (DIS) 

8.1. Bjorken sum rule 


7.2. Higgs decay into quarks 

The (inclusive) decay width of the Higgs boson into 
a pair of quarks is related to the spectral density of the 
scalar correlator according to (3|5 1 . For ny - 5 which 
corresponds to the newly discovered Higgs boson we 


The Bjorken sum rule (for polarized DIS) expresses 
the integral over the spin distributions of quarks inside 


^We have extended the package by including the five-loop effects 
to the running of Uj and quark masses. 
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Figure 3: Examples of diagrams contributing to the coefficient func- Figure 4: Examples of diagrams contributing to the coefficient func¬ 
tion at three and four loops. tion at three and four loops. 


of the nucleon in terms of its axial charge times a coef¬ 
ficient function 


rr(e') 



Q^) - gf (.r, Q^)]dx 


o 


Qli-l ’ 
i=l ^ 


(74) 


where g"y and g™ are the spin-dependent proton and 
neutron structure functions, is the nucleon axial 
charge as measured in neutron y6-decay. The coefficient 
function C^^P{as) - 1 +0{as) is fixed by the OPE of two 
EM currents (for a more detailed discussion, see If84l l: 


f 


T[J^(x)J^(0)]e‘^’^dxl 


(75) 


9 


,2 ■'■>,) A;( 0 ) + ... 


Ca^'^iOs) - (76) 

Tr[E\] + Tr(£) Tr[£fJ Cff (a,) = 

Jc^f(a,) + 3Tr[£] C^f (fl.)j Tr[E\] 

= C‘^jP(a,)TT[E\]. 


Here E = diagiQi) is the quark charge matrix, = 
il/Ejpili is the quark EM current, - tf/jp if/ is the 
(flavor non-singlet!) axial current and - —q^. 

As one can see from eqs. (|76]l, the coefficient function 
C’^jP{a,) receives contributions from two types of dia¬ 
grams, viz. the singlet and non-singlet ones (see Figs. 
and [^respectively). Let us discuss them in turn. 




+ a. 
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24 

103 ^ 1 2 5 X 

+ ^^3 - ^ ^3 + 7 ^ 6 ) 


+ n 
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f 5832 


(77) 


= 1 - ( _ 4.583 -h 0.3333 n/) (78) 

-I- a3 ^ - 41.44 -H 7.607 «/ - 0.1775 nJ) 

-b a;J(-479.4-H 123.4«/-7.697 n^H-0.1037 n}). 

Note that for phenomenologically relevant values of 
<3 GeV^ one should work in an effective QCD with 
three active flavours. In this case the singlet contrubu- 
tion vanishes identically as Tr[£’] s 0 and 

C'^jPinf ^3) = C%P(nf = 3) = (79) 

I-a,- 3.583^3 - 20.22al - 175.7 


The coeffcient function C^-'Pias) starts from 1 which 
corresponds to the parton approximation. The four-loop 
result was published in ll85l for the case of a general 
gauge group. The QCD result reads: 


Phenomelogical implications of (79 1 have been studied 
in works ll86l l87l . Their results can be summarized as 
follows. 
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Figure 5: Perturbative part of the Bjorken sum rule as a function 
of the momentum transfer squared in different orders against the 
combined set of the Jefferson Lab (taken from V.L. Khandramai, R.S. 
Pasechnik, D.V. Shirkov, O.P. Solovtsova, O.V. Teryaev, Four-loop 
QCD analysis of the Bjorken sum rule V5 data, Phys.Lett.B706:340- 
344,2012). 


First, by comparing experimental data (see Fig. 
with the theoretical prediction ( [79| the authors have ar¬ 
rived at the following conclusion ll8^ : ’’One can see 
that at > 0.7 GeV^ the four-loop approximation de¬ 
scribes the data quite well. Moreover, the corresponding 
curve passes close to the central values of several data 
points, although the experimental accuracy (which is of 
the same order as both the three- and four-loop contri¬ 
butions) does not allow one to make a definite choice 
between four- and three-loop approximations.” 


Second, a certain duality between higher orders and 
higher-twist contributions to the Bjorken sum rule (74i 
has been detectecQin ill. Indeed, the fitted value of 
first non-zero higher twist contribution, yU 4 , has proved 
to be strongly dependent on the order of PT terms kept 
in For example, at leading order (that is with all 
terms in except first two set to zero) //4 was found 
to be -0.037 + 0.003 GeV^. At next-leading-order /i 4 is 
decreased to -0.025+0.04 GeV^ and, finally, at (9(a^) it 
becomes compatible to zero; — 0.005 + 0.008 GeV^. 


The singlet contributions to formally starts at 
two loops but the corresponding diagrams sums to iden¬ 
tical zero due to Furry’s Theorem. The next, three-loop 
term also happen to be zero ||2T1 . This fact has been 
explained in Il84ll with the help of generalized Crewther 
relation ll89H^ . Here it has been also predicted that at 
four loops the singlet contribution should have the form; 


^The phenomenon in a more general context was earlier discussed 

in 1551 . 





(80) 


withySo = \jCa - = 2Tr({f, y)f)), and X 

being a constant. 

We have performed a direct calculation of C^f at or¬ 
der aj. Our result reads; 


in full agreemenj^with (801. 


9. Anomalous Dimensions of twist 2 operators 

Recently there has been a lot of progress in three-loop 
QCD calculations of the moments (and the correspond¬ 
ing anomalous dimensions) of deep inelastic structure 
functions Il92lj96l ). In particular, the anomalous dimen¬ 
sion of of the twist-two non-singlet operator (i// 
and ijj' refers to the two different quark species) 

has been analytically found for a generic value of spin 
N at the three loop level. 

In fact, a general consistency argument requires the 
use of four-loop splitting functions in applications of 
the results of ll94l l96l to the phenomenological analy¬ 
sis of deep inelastic experimental data. Unfortunately, 
the formidable problem of tht four-loop calculation for 
generic N seems to be out of reach for available tech¬ 
nologies. On the other hand, fixed N calculations are 
now possible (at least for not too large values of N). 

The first result in this direction was reported in ll97l 
where the four loop anomalous dimension of the oper¬ 
ator was computed for a particular number of 

quark species ny — 3. Later this result was confirmed 
and generalized to a generic gauge group in ||98l. Note 
that at four-loop level the calculation of y’^^ with the 
use of BAICER does not require application of any IRR; 
one just computes a diagonal matrix element 

with Ip) being an off-shell quark state. In principle one 
could even compute the five-loop anomalous dimension 


^In fact, paper EH has also guessed a particular value of 
X = - 5 ( - ^ + ^^3 “ ^^5 ) which happens to be very different 
from our result. 
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for low N. However, this would require signifi¬ 
cantly more computer as well as human power (the lat¬ 
ter due to quite complicated IRR). 

Two years ago the present authors computed for 
N - 2,3 and 4 in QCD with full tip dependenc^ Our 
results read; 
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^ The results given below were first presented on the 19 Meeting 
of SFB/TR9 “Computational Particle Physics” 19.03.2013 (Aachen). 
Very recently and have been computed for a case of a generic 
gauge group 1991 . For the QCD case gauge we have found full agree¬ 
ment between eqs. (82]93) and results of ED- 
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Numerically all 3 anomalous dimensions display a re¬ 
markable similarity (modulo global normalization): 

y^Hrif = 3) = ^ (a, 4- 2.73190^ -H 7.8763a^ 4- 28.706a^) , 

y^Hrif = 3) = ^ (a, 4 - 2.49820? -H 7.0891a? 4 - 23.587a?) , 
yiHrif = 3) = (a, + 2.3871a? -H 6.8288a? 4- 22.294a?) . 
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10. Conclusions 

The problem of analytical evaluation of massless 
propagators at four loops has been under investigation 
since long Ilia. It has been solved using reduction via 
l/D expansion. As a result a number of important four 
and five loop calculations have been done. In this short 
review we have briefly discussed some of them related 
to the R-ratio, Higgs decays into quarks, deep inelastic 
scattering and QCD renormalization group functions. 
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